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We calculate all multipoint correlation functions of all local bond modifications in the two- 
dimensional Abelian sandpile model, both at the critical point, and in the model with dissipation. 
The set of local bond modifications includes, as the most physically interesting case, all weakly 
allowed cluster variables. The correlation functions show that all local bond modifications have 
scaling dimension two, and can be written as linear combinations of operators in the central charge 
—2 logarithmic conformal field theory, in agreement with a form conjectured earlier by Mahieu and 
Ruelle in Phys. Rev. E 64, 066f30 (2001). We find closed form expressions for the coefficients 
of the operators, and describe methods that allow their rapid calculation. We determine the fields 
associated with adding or removing bonds, both in the bulk, and along open and closed boundaries; 
some bond defects have scaling dimension two, while others have scaling dimension four. We also 
determine the corrections to bulk probabilities for local bond modifications near open and closed 
boundaries. 



PACS numbers: 05.65.+b,45.70.-n 



I. INTRODUCTION 

Self-organized criticality may be the underlying cause 
of power laws in a wide range of natural and man-made 
phenomena [1, 2]. Systems exhibiting self-organized crit- 
icality naturally approach a critical state, without any 
intrinsic time or distance scales. The critical point is 
reached without any fine-tuning of parameters. This is 
unlike most critical points seen in physics — for example, 
the Ising model is only critical at a single, very specific, 
temperature. This lack of fine-tuning is essential if we are 
to understand power laws in nature, where no fine-tuning 
is possible. 

Since the concept of self-organized criticality was first 
introduced by Bak, Tang, and Wicscnfcld in 1987, a num- 
ber of models have been developed to investigate this 
phenomenon [3]. However, the original model, the two- 
dimensional, isotropic, Abelian sandpile model (ASM), 
is still one of the simplest and most interesting of the 
models. The ASM is simple and robust, which are nec- 
essary features for any model of self-organized criticality. 
While natural phenomena are quite complex, any model 
that seeks to explain the ubiquity of power laws in na- 
ture must, paradoxically, be very simple; if we are to 
have a robust model for the generation of power laws, we 
must neither have finely-tuned parameters, nor finely- 
tuned rules. 

The ASM is defined on a lattice of sites, and is de- 
scribed by a toppling matrix A, whose dimension is equal 
to the number of sites in the sandpile. The sandpile 
evolves stochastically. In each time step, a grain of sand 
is added to a random site. Then, sites are checked for sta- 
bility. If the number of grains at a site i is greater than 
A^ > 0, then the site i is unstable, and topples, losing 

A? y grains, while every other site j gains -Ay j > grains. 
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(Generally, Ay j is zero except when j neighbors i.) Typ- 
ically, models are conservative, which means that each 
toppling in the bulk conserves the total number of grains 
(^jAjj = 0). Only for topplings along the boundary, 
where grains can fall off the edge, can the total number of 
grains change. We continue toppling unstable sites until 
no sites are unstable. Then, wc begin a new time step, 
and again add a grain to a random site. 

The ASM is surprisingly tractable [4-6]. We only 
briefly cover some of the essential points here — for com- 
prehensive reviews, see [7, 8]. 

After a large number of time steps, the ASM reaches 
a well-defined distribution of states. Of the stable height 
configurations, some are transient, and occur with prob- 
ability zero after a long amount of time. All other states 
are recurrent, and occur with equal probability. Dhar 
showed that the total number of recurrent states is just 
det(A) [5]. This is also equal to the number of spanning 
trees that can be drawn on the lattice, showing a connec- 
tion between the sandpile and spanning tree problems [6] . 

These statements hold for all ASMs, which define a 
large class of models. Now, we specialize to the two- 
dimensional, conservative, isotropic ASM, which is de- 
fined on a two-dimensional square lattice, where each site 
has a maximum height of four, and where upon toppling 
at any site, one grain is sent to each of the site's four 
neighbors. Furthermore, we work in the limit where the 
lattice is infinite. The two-dimensional, isotropic, span- 
ning tree problem is equivalent to the central charge -2 
logarithmic conformal field theory (c = —2 LCFT) [6], 
which has the simple Gaussian action S = (1/V) J 8989, 
where 9 and 9 are complex Grassman variables. The 
c = -2 LCFT is described in [9-11]. While the two- 
dimensional, conservative, isotropic ASM is just one 
of many possible ASMs, it is the original, standard, 
model [3], and it is reasonable to simply refer to it as 
"the ASM," which we do for the remainder of this paper. 

Calculations of correlation functions, using methods 
to be described in the next section, have confirmed that 
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there is a relationship between the ASM and c = — 2 
LCFT. Two-point correlation functions of unit height 
variables decay as 1/r 4 in the bulk [4], as do all two- 
point height correlations along open and closed bound- 
aries [12]; these correlations can be understood as equiv- 
alent to correlations of LCFT operators. Furthermore, 
calculations of certain three-point correlation functions 
of heights along closed boundaries, and all multipoint 
correlations of heights along open boundaries, have al- 
lowed LCFT field identifications for heights along bound- 
aries [13, 14]. 

The ASM is not robust to all perturbations. If wc 
relax the constraint that the model be conservative, and 
instead allow grains to be lost in any bulk toppling (i.e. 
allow dissipation), correlations decay exponentially, and 
we are taken off the critical point [15-17]. The condition 
of conservation can be considered a "natural" one, rather 
than one requiring "fine-tuning." Deeper probes of the 
the conformal structure can be obtained by looking at 
correlations both on and off the critical point. 

Mahieu and Ruelle calculated a number of off-critical 
correlation functions of certain height configurations, 
known as weakly allowed clusters (WACs) , and used their 
correlation functions to propose field identifications for 
the 14 simplest WACs [18]. They found that their corre- 
lation functions could be explained by assuming that all 
14 WACs took the form 

<j){z) = -{a-, oede + deae ■. +b 1 ■. oeoe + dede ■. + 

+iB 2 : 8986-8959 : . qq A (1) 

2ir J 

The coefficients A, B 1 , B 2 , C, and P(S) vary from WAC 
to WAC. P(S) is the probability for the cluster at the 
critical point, and M is the mass, a measure of how 
far the model is from the critical point. The correla- 
tion functions that they used were mostly two-point func- 
tions along horizontal or diagonal well as some 
three-point and four-point functions for the two simplest 
WACs. 

While these calculations provide strong evidence for 
the identification of the ASM with the c = -2 LCFT, and 
the field identification in Eq. (1), the fact that only spe- 
cific correlations were considered limits the range of the 
identification. It would be surprising if new orientations 
of correlation functions, or new WACs, were found to 
be inconsistent with Eq. (1); but the calculations in [18] 
do not rule this possibility out. More importantly, since 
each correlation function in [18] required a new and sep- 
arate calculation, it is hard to understand, mathemati- 
cally, why these results occured. While their end results 
showed that certain correlations of WACs in the ASM 
are equal to correlations of Eq. (1) in the LCFT, it was 
not mathematically transparent as to why this should be. 
Nor was it clear why, or if, the same coefficients would 
appear in other properties, such as off-boundary correla- 
tions, or correlations with defects. 



Here, we calculate all correlation functions of all lo- 
cal bond modifications (LBMs), for arbitrary numbers 
and types of LBMs at arbitrary positions (far from one 
another); our calculations confirm that all LBMs should 
receive the field identification in Eq. (1). By LBMs, we 
mean any set of local changes in the sandpile toppling 
rules. For the ASM at the critical point, we will as- 
sume that the LBMs are conservative (do not create or 
destroy grains). Since all WACs can be calculated by 
LBMs, our results automatically include all correlations 
of WACs. While the WACs are the most important types 
of LBMs, and the easiest to find probabilities of, in nu- 
merical simulations, we generally discuss our results in 
terms of LBMs, to emphasize the generality of our re- 
sults. We give closed form expressions for A, B\, B 2 , 
and C, and describe methods that allow rapid calcula- 
tion of these coefficients. While a computer is needed for 
the calculation of specific A, B\, and B 2 coefficients, the 
general calculations can be done by hand. 

By showing how calculation for all LBMs can be done 
at once, we make the mathematical structure clearer. For 
example, we can quickly see why the coefficients A, B\, 
and B 2 appear in other properties. We illustrate this by 
looking at off-boundary LBM probabilities, and correla- 
tions with bond defects (either in the bulk, or along a 
boundary). Interestingly, we find that some bond defects 
are represented by a LCFT operator with a scaling di- 
mension of four; this is, to the best of our knowledge, the 
first dimension four operator found in the ASM. 

While our calculations have been done in both the nor- 
mal ASM, and the ASM with dissipation, wc focus our 
discussion on the simpler analysis at the critical point, 
and only discuss the more complicated massive correla- 
tions in the last section, and in the appendices. 



II. WEAKLY ALLOWED CLUSTER 
VARIABLES 

The methods used in this paper are not powerful 
enough to calculate probabilities and correlations for any 
height configurations. Even the calculation of the prob- 
ability for a site to have height two requires much more 
complicated methods [19], and the correlation function of 
two height two variables remains unknown. This is be- 
cause the condition for a site to have height two involves 
a nonlocal condition. 

As already stated, the most important LBMs are those 
used to calculate properties of WACs [20] . WACs are re- 
lated to forbidden subconfigurations (FSCs). An FSC 
is a height configuration over a subset of sites F, such 
that for every site i € F, the number of neighbors of i in 
F is greater than or equal to the height at i. FSCs are 
important because ASM height configurations are recur- 
rent if and only if they have no FSCs [5]. A WAC is a 
height configuration that contains no FSCs, but becomes 
a FSC if any height in the WAC is decreased by one. 
Three WACs are shown in the left side of figure 1 . 
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FIG. 1: Some WACs and their corresponding modified sand- 
piles 



WACs are analytically tractable because it turns out 
that the number of sandpiles with a particular WAC is 
equal to the number of recurrent states in a sandpile with 
modified toppling rules [4]. There are actually several 
different ways to modify the toppling rules to obtain the 
WAC probability. The simplest is to, for each connected 
piece of the WAC, remove all but one of the bonds con- 
necting it to the rest of the lattice — the modified lattices 
corresponding to the WACs are shown in the right side 
of figure 1. (See [18] for a discussion of other ways in 
which the sandpile can be modified to obtain the WAC 
probabilities.) In these modified sandpiles, grains of sand 
cannot flow along the removed bonds; to continue to con- 
serve the number of grains during each toppling, the con- 
dition for instability must be decreased at the sites at the 
end of the removed bond. These changes result in a new 
toppling matrix A'. 

As already stated, the number of recurrent states in 
the ASM is det(A). The number of recurrent states that 
have the WAC is given by det(A'). (We discuss this 
equivalence further in section V.) So the bulk probability 
for the WAC is given by 



P 



dct(A') 
det(A) 



det(I + BG) , 



(2) 



where we have defined B = A' — A and G = A -1 . G 
is the well-studied lattice Green function (at the critical 
point); exact expressions are known for the Green func- 
tion between nearby sites, and asymptotic expressions 
for the Green function between distant sites [21]. While 
A, A', and G all have large dimensions (equal to the 
number of sites), B is zero outside of a finite collection of 
sites. When the bond between i and j is removed, B^j 
and B? ? arc both increased by 1, while B^r and B? 7 are 
both decreased by 1. For example, for the unit height 
probability, we have 



32 33 



B 




(3) 



Here i is the site fixed at height one, while ji, j'2, and j'3 
are the three sites that i has been disconnected from. 

For any WAC, the fact that B is finite-dimensional 
means that the height probability can be found by calcu- 
lating a simple, finite-dimensional, matrix determinant. 
All WACs thus correspond to LBMs. However, many 
LBMs do not correspond to WACs. LBMs are simply any 
sandpile modifications that can be modeled with a B ma- 
trix that is conservative (every row and column sums to 
zero) and symmetric. Our analysis gives all correlations 
of LBMs, which thus automatically gives all correlations 
of WACs. 



III. CORRELATIONS OF LOCAL BOND 
MODIFICATIONS 

For an n-point correlation function of LBMs, we can 
still use this method. The only difference is that the 
removed bonds are located in n distant clusters; this is 
illustrated in figure 2. Removal of bonds in this fashion 
will give B and G block matrix structures. For example, 
for a 3-point function, we will have 







\ 







B 2 




V 


B 3 J 




(Gil 


G12 G13 




G2I 


G22 G23 
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G32 G33 



(4) 
(5) 



B M is the modification to the toppling matrix for the set 
of bonds removed about the u th LBM. G uu is the Green 
function matrix between sites of the w th LBM, and its 
elements are 0(1). G uv , u ^ v is the Green function 
matrix between the sites of the u th and v th LBMs, and 
its elements are given by the bulk Green function, 



G (x,y) = -—ln(x 2 +y 2 ) 



J_ 

2tt 



In 8 
4ir 



+ ... 



(6) 



where 7 = 0.57721 ... is the Euler-Mascheroni con- 
stant [21]. We work in the limit where the LBMs are all 
very far from each other — we assume that any two of the 
n LBMs are the same order-of-magnitude, O(r), apart. 
Since the Green function diverges as ln(r) with increasing 
r, calculation of det(I + BG) initially looks very difficult. 
However, every row of every B u sums to zero — this fol- 
lows from the manner in which we constructed B M , and 
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(r uv cos 4> uv , r uv sin <j> uv ) . The w th LBM covers a set of 
sites at locations (ki,h), relative to (0,0), and the V th 
LBM consists of a series of sites at locations {k 2 , h), rela- 
tive to (x uv ,y uv ). (ki, k2, h, and l 2 are all 0{1).) Then, 
the elements of G uv all have the form Go(x uv + k 2 — 
k\,y U v + h — h)- The last two paragraphs show that 
we only need the parts of G uv that depend on both the 
row and column indices. That is, we only need the parts 
of the Green function that depend on both (ki,h) and 
(k 2 , h), and can drop all other terms. Expanding Eq. (6) 
in powers of l/r uv , we find that the lowest-order term 
not dropped is 



FIG. 2: Modified sandpile for a WAC two-point correlation 

reflects the fact that grains of sand are still conserved in 
each toppling in the bulk of the modified sandpile. This 
implies that parts of G uv that depend only on the column 
index make no contribution to BG, and thus no contri- 
bution to the correlation function, det(I + BG). So we 
only care about differences (discrete derivatives) of Green 
functions between columns of G uv , and the elements of 
G uv are effectively 0(l/r), rather than O(lnr). 

For LBMs, every B u is symmetric, so every column 
of every B M sums to zero. Using the matrix identity 
det(I + BG) = det(I + GB), this in turn means that the 
parts of G uv that depend only on the row index make 
no contribution to the probability. This is, in effect, like 
taking another discrete derivative of the Green function, 
so that the elements of G uv are effectively 0(l/r 2 ). 

To make this concrete, suppose that the local ori- 
gin of the u th LBM is located at (0,0), and the lo- 
cal origin of the v th LBM is located at (x uv ,y uv ) = 



G (x uv + k 2 - fci, y U v + h - h) -» --z — — x 

^ 7rr uv 

((fcifc 2 - hh) cos(2<j) uv ) + (fciZ 2 + k 2 k) sin(20 ut) )) 

(7) 

The Green function can thus be treated as 0(l/r 2 ) for 
correlations of LBMs. For more general local arrow di- 
agrams, such as those that appear in the calculations 
involving the height two variable [12, 19], the B matrices 
are not symmetric, and we can no longer drop the parts 
of G uv that depend only on the row index. 

To get the connected n-point function from det(I + 
BG), we need to pick at least one element off the block 
diagonal in every block row of G, and in every block 
column of G. This means, at the minimum, picking 
n elements off the block diagonal of G, resulting in a 
leading-order contribution to the correlation function of 
(D(l/r 2n ) — this is the universal part of the correlation 
function. 

Mahieu and Ruelle showed that for two-point func- 
tions, the constraint of picking only two elements off the 
block diagonal allows the correlation function to be writ- 
ten as [18] 



det(I + BG) = -p ul p U2 Trace I B Ml G MlM B„ 2 G„ 2Ul \ (8) 

Similarly, they found that the leading-order contribution to the three-point probability is 



det(I + BG) = 

f I I II 

PuiPu 2 Pu 3 Trace < _ — T5 Ul G Ul u 2 ^ R 7^ B U2 G M2 „ 3 — B U3 G„ 3Ml > + 

PuiPu 2 Pu 3 Trace < ,-r> 7^ B Ul G UlU3 _ _ B U3 G M3M2 pq B U2 G„ 2Ml > (9) 



More generally, for an n-point correlation, if only n terms are picked off the block diagonal, then the connected 
correlation function is given by 



■5 



det(I + BG) 



(-) 



n+1 



,x=l 



^ Trace I] 



1 + B„ G 



(as) "".(^"sf.) 



(10) 



where s is summed over all one-to-one mappings from 
{1, 2, ... n} to {1, 2, . . . n}, and we identify s(n + 1) with 
s(l). Mahieu and Ruelle wrote Eqs. (8-9) in different, 
but equivalent, forms. 



p u = dct(I + B U G UU ) 



(11) 



is the bulk probability of the u th LBM. (Note that the 
two trace terms of Eq. (9) are actually equal. We have 
written the three-point function in this form to make 
clear how the form generalizes for n-point functions.) 
We can rewrite Eq. (7) as 



Cu,kk — Puk u 
Cu,kl 



I 



-PuZ 



I 

1 I 



c u ,n = -pj 



rr 



I + B U G U „ 



-B,, 



(15) 

(16) 
(17) 



Then, the correlation function of n LBMs is given by 



G„ 



27T7-2 



{{kukl - coa(2cj) uv )- 



(WT + C^)sin(2^)) (12) 



k u is the column vector of the horizontal positions of the 
sites of the u th LBM, relative to the u th local origin (i.e., 
the elements of k u are 0(1)). l u is the corresponding 
vector of vertical positions. k u and l u are both length 
N u , where N u is the number of sites needed to represent 
the u th LBM with the methods of the previous section 
(e.g. N u = 4 for the unit height variable). 

We insert Eq. (12) into Eq. (10). For each of the n 
G ra 's, we can pick any of the four matrices of Eq. (12), 
resulting in 4™ terms. In each of these 4" terms, each G uv 
has been replaced with the product of a column vector 
and a row vector. Using the cyclicity of the trace to 
move one row vector at the end of the trace to the start 
of the trace, we see that each matrix (I+B„G tlu )~ 1 B tt is 
bracketed by a row vector to its left, and a column vector 
to its right, producing a 1 x 1 matrix. So each of the 4™ 
terms is the product of n numbers. We can represent 
the decisions as to which terms of Eq. (12) to pick by 
representing G uv with a 2 x 2 matrix, N Mt) . The possible 
ways to bracket (I + B u G uu ) _1 B tl can be represented 
with a 2 x 2 matrix, M u . We have 



k„ 



M„ = 



N = 



Cu,kk C-u,kl 



(13) 



rr 



1 



We have defined 



cos(2(j) uv ) sin(20 ut) ) 
sm(2c/) uv ) -cos(2<p uv ) 



k u 



(14) 



-Trace (M U1 N„ 1U2 M„ 2 N 

U n Ul ) 

— (((n — 1)! — 1) other trace terms) , (18) 

where the other trace terms are derived by permutations 
of {ui, U2, ■ ■ ■ , u n }, as in Eq. (10). 

We can compare this to correlation functions of fields 
in the c = — 2 LCFT. Mahieu and Ruelle proposed that 
the WACs are represented, at the critical point, by 



<jy u (z u ) = -{A u : 8689 + 8989 : +B lu : 8989 + 8989 : 
+iB 2u : 8989 - 8989 :} (19) 

(The "C" term in Eq. (1) only appears off the critical 
point.) 

We can compute connected n-point correlations of 
these fields in the c = — 2 LCFT. We use the formula- 
tion of the c = — 2 LCFT where the action is 



8989 : 



(20) 



where we don't integrate over zero modes in expectation 
values. Since the theory is Gaussian, to calculate correla- 
tion functions we simply need to take Wick contractions. 
The relevant nonzero ones are 



(89(z u )89(z v )) 
(89(z u )89(z v )) 



-2id 



l 

2(z u z v )^ 



2r 2 

,+H4>uv 



2r 2 



(21) 



(22) 



Each term of Eq. (19) has one 9, and one 9. The 
only difference between terms is whether the derivative 
on the 9 is holomorphic or antiholomorhic, and whether 
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the derivative on the is holomorphic or antiholomor- 
phic. We can use a 2 x 2 matrix to represent the choice 
of which terms of <j) u {z u ) are picked: 

de Be 

„ de ( B lu + iB 2u A u \ 
u ~ Be { A u B lu -iB 2u ) [Z6} 

The contractions of Eqs. (21-22) can then be represented 
with the matrix 

de Be 

de (- e-™™/(2rl v ) \ 

UV ~ &e\ -e+ 2 *W(2rl) ) 

(24) 

The contribution to the correlation function where the # 
of the 1st LBM contracts with the 6 of the 2nd LBM, 
the e from the 2nd LBM contracts with the of the 3rd 
LBM, and so on, is 



calculation can be made substantially faster with the fol- 
lowing matrix identity, which we state without proof: 

det (l + B (G + Kff^ = dct(I + BG) + 

+Kdet (I + BG) (> f -L_B/) (29) 

This identity holds for any vectors / and g, and c- number 
K. It allows us to compute the c's in Eqs. (15-17) as ma- 
trix determinants, which is faster than computing matrix 
inverses. Furthermore, we note that, in general, certain 
combinations of rows (and columns) of B M will sum to 
zero, which means that we can perform a matrix rotation 
to reduce the size of the matrix determinant. With these 
methods, computation of A, B\, and B 2 for the ten sim- 
plest WACs takes roughly one hour, using Mathematica 
on a computer with a 1.2 GHz processor. The results 
agree with those found in [18]. Comparison with two of 
larger WACs, which they label Sio and Sn, requires a 
more detailed discussion of the mapping between WACs 
and LBMs, which is done in the next section. 



Trace(F Ul H t(1 U2 F U2 H U2U3 . . . F Un H UnUl ) (25) 

Other contractions give other permutations, just as in 
Eq. (18). Finally, M„N„„ differs from F U G UV only by 
a matrix rotation, which will not affect the trace, if we 
make the following identifications: 

A u = — (c u .kk + Cull) (26) 

Bin = -z-(c u .kk — Cu.u) (27) 

B 2u = —Cu.u (28) 

So the traces in Eq. (18) and Eq. (25) are equal, and 
all LBMs are indeed represented by the field in Eq. (19). 
These formulas for the coefficients have the appropriate 
transformation properties under 90° rotations, and x and 
y reflections. (Technically, the overall sign of Eq. (26) 
is still undetermined at this point, since all correlation 
functions have even numbers of A's. To determine the 
signs of the A's we need to look at at least one massive 
correlation function. We can do this by consulting the 
massive three-point function of the unit height variable 
in [18], or more broadly, by looking at the general massive 
correlations in section IX.) 

IV. COMPUTATION OF A, B x AND B 2 TERMS 

A, B\, and B 2 , can be calculated on a computer with 
Eqs. (26-28) and Eqs. (15-17). Evaluating Eqs. (15-17), 
as written, requires taking a matrix inverse, which can be 
computationally time-consuming for larger LBMs. The 



V. THE MAPPING BETWEEN WEAKLY 
ALLOWED CLUSTER AND LOCAL BOND 
MODIFICATIONS 

We illustrate the mapping between larger weakly al- 
lowed clusters and local bond modifications with the 
sandpile modification shown in figure 3. In this modi- 
fied sandpile, a five-site cluster is separated from the rest 
of the sandpile, except by a single bond. The number of 
states in the modified sandpile of figure 3 is equal to the 
number of states of the unmodified ASM where decreas- 
ing the five-site cluster's leftmost site (which we call T) 
from 2 to 1 makes the five-site cluster a FSC, and does 
not produce any larger FSCs [19]. The condition that the 
FSC produced be maximal is necessary for this equiva- 
lence, although this condition was not explicitly stated 
in [19]. If changing the height of i from 2 to 1 makes 
the five-site cluster a FSC, the original height configura- 
tion (before this change) must have been one of the four 
configurations shown in figures 4 and 5. Of these four 
configurations, the one in figure 4 is not a WAC, while 
the three in figure 5 are. 

Absent other conditions, the configuration in figure 4 
does not have the same probability as the configurations 
in figure 5. The configuration in figure 4 is more likely, as 
an allowed configuration always stays allowed when a site 
height is increased. However, for some configurations of 
heights outside the five-site cluster, decreasing one of the 
three height two sites other than the one at i, to height 
one, will create an FSC outside the five-site cluster, so 
that figure 4 is allowed, but figure 5 is not. 

However, if we impose the condition that taking the 
height of i from two to one should produce no FSC larger 
than the five-site cluster, then all four configurations in 
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4^ 



t r r 



FIG. 3: Modified sandpile for the height configurations in 
figures 4 and 5 



<2H2)- 



FIG. 4: A non-WAC with the same shape as the configuration 
in figure 3 



with theirs. 

Although wc chose a specific sandpile modification, the 
discussion is easily generalized. Generally, consider sand- 
pile modifications similar to those in figure 3, which sepa- 
rate a cluster of sites from the rest of the sandpile, except 
for one linking site. There will be N possible height con- 
figurations in the cluster that become FSCs if the height 
of the linking site is reduced. Some of these will be WACs 
and some will not. The probability associated with the 
lattice modification will be N times the probability for 
any individual WAC. (It will also be times the prob- 
ability for any of the other height configurations, given 
the condition that decreasing the height of the link site 
should produce no FSC larger than the cluster in ques- 
tion.) 

These issues did not need to be discussed for the sim- 
pler WACs, such as those in figure 1 (or the other WACs 
of [18]). For each of the sandpile modifications in the 
right side of figure 1 , only one corresponding WAC height 
configuration is possible (N = 1). 



figure 4 and 5 are equally probable. With this condition, 
decreasing one of the height two sites in figure 4 cannot 
possibly create an FSC outside the five-site cluster, since 
then the union of this FSC with the five sites would be a 
larger FSC upon decreasing the height of i. 

For a WAC, it can be shown that if any height is de- 
creased, not only does the WAC become a FSC, but it is 
not contained in any larger FSC. Therefore, for the three 
height configurations in figure 5, the condition that the 
five-site FSC generated is maximal is automatic. (How- 
ever, for the configuration in figure 4 it is not.) Therefore, 
the probability associated with the modified sandpile in 
figure 3 is four times the probability of any of the three 
WACs in figure 5. It is also four times the probability of 
the height configuration in figure 4, if we impose on this 
configuration the condition that decreasing the height of 
i should produce no FSC larger than the five-site cluster 
(although clearly this is not as physically interesting). 

In [18], the configurations that they labelled Siq and 
Sn were not WACs. Once they are modified to be WACs 
with the same shape, our values for A, B\, and B 2 , ob- 
tained with the methods of the previous section, agree 




-(2X2)- 



I 




-(2)^3)^1 



FIG. 5: Three WACs with the same shape as the configuration 
in figure 3 



VI. OFF-BOUNDARY EXPECTATION VALUES 

These methods allow us to quickly determine the ef- 
fects of a number of defects or boundary conditions on all 
LBMs, and see that each time, we obtain the same coef- 
ficients. In many cases, a defect or boundary changes the 
Green function matrix in a manner such that the change 
factorizes, taking the form Kfcf, for some vectors / and 
g. Eq. (29) then shows that the effects of the change can 
be written as a linear combination of the c's, and thus as 
a linear combination of A, Bi, and B 2 . 

For example, consider the probability for an LBM lo- 
cated a distance y from a boundary (open or closed). 
Let the boundary be at y — 0, and x be the coordinate 
along the boundary. Then the Green function is modified 
to [22] 



G open {xi,yi;x2,yi) = G (xi - x 2 , 2/1 - yi) 

-G (zi - a*,!/! +y 2 + 2) (30) 
G c i OS ed (£1,2/1; x 2 ,y 2 ) = G (xi - x 2 ,yi - 2/2) 

+G (xi - x 2 , 2/1 + 2/2 + 1) (31) 

Placing the local origin of the LBM at (0,y), the Green 
function between points (fci,Zi) and (k 2 ,l 2 ), relative to 
this local origin, where fci, h, k 2 and l 2 are all 0(1), is 



G(ki,h;k 2 ,l 2 ) = G {ki - k 2 ,h - l 2 ) ± 



+0 



k\k 2 + l\l 2 



8ny 2 



y- 



(32) 



In the ±, the top sign is for closed boundaries, and the 
bottom sign is for open boundaries. As with Eq. (7), we 
have only kept terms that depend on both (k\,li) and 
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(k 2 ,l 2 ) (otherwise, there would be terms of 0(lny) and 
(D(l/y)). Then, using Eq. (29), keeping only terms of 
0(l/y 2 ), and using the definitions in Eqs. (26-28) and 
Eqs. (15-17), we immediately see that the probability for 
the LBM is 



This is consistent with identifying the bond defect with 



fcbond 

2tt 



: (q x d x 6 + q y d y 9)(q x d x e + q y d y 6) : (39) 



y3 



(33) 



This agrees with results for the unit height variable found 
in [22]. 



VII. BOND DEFECTS IN THE BULK 



VIII. BOUNDARY OPERATORS AND BOND 
DEFECTS 

It was shown in [14] that any local arrow diagram along 
an open boundary is represented in the LCFT by the 
operator 



We have also used the methods described here to inves- 
tigate bond defects in the ASM. We change the toppling 
matrix from the defect-free matrix Aq to 



A(m) 



SA(i,j) 



A {i,j) + 5A(i,j) 

-fcbond if = (so, so) or 

(hi) = (si, si) 
+fcbond if 0J)j= (so, si) or 
(*, j) = (si, s ) 



(34) 



(35) 



If so and si are adjacent sites, and fcbond = 1, then this 
corresponds to removing the bond between s and si . If 
fcbond = — 1, this corresponds to adding a bond between 
s*o and si. 

The Green function is the inverse of the toppling ma- 
trix, and the effects of this perturbation can be calculated 
by summing a geometric series. The result is 



■- det(I + B U G UU ) ((y + Tf—L—B u (y + 1) 



xd6d9 , 



(40) 



where y is the vector of distances perpendicular to the 
boundary. The arguments there worked for any local 
arrow diagram (not just those corresponding to LBMs), 
and similarly to our arguments in section III, used the 
fact that the Green function along open boundaries falls 
off as 1/x 2 . Along closed boundaries (and in the bulk), 
the Green function grows as lnx, so the situation is more 
complicated, and not all operators are proportional to 
8889. (See, for example, the operators for the height 
two and three variables along closed boundaries, given 
in [13, 14].) However, we can now derive an expression 
similar to Eq. (40) for LBMs along closed boundaries. 

The Green function for two sites on a closed boundary, 
and x » 1 apart, was found in [22], using Eq. (31), to be 



G(hj) = G (lJ) + fcbond (G (?, s ) - G (i, si)) x 

(Go(s-*o,J)-Go(si,J)) ,(36) 

where 

1 1 

fcbond fcbond 

G (s , s ) + Go(si,si) - G (s , si) - G (si, s ) 

(37) 

The correction to the Green function factorizes, just as 
in Eq. (29), and the corrections to the LBM probabilities 
again depend on the same coefficients. If the two ends of 
the bond defect are at (0,0) and (q x ,q y ), where q x and 
q y are both 0(1), then from Eqs. (15-17), Eqs. (26-28), 
and Eq. (29), the probability for a LBM of type u at 
(x,y) — (r cos 6*,r sin 9), for r » 1, is 

fcbond r / 2 , 2 \ a 

~4^t {(q x + q v )A u + 

(ql - q 2 y)(Bm cos(49)+B 2u sm(49j) + 
{2q x q v ){B lu sin(40) - B 2u cos(40))} (38) 



1 , , (l ln2\ 1 
G do , c „(0,0;,.0) = --,„(,)- (2 + _j+_ 

< 41 > 

Using the recursion relation G A = I on this equation, we 
can extend Eq. (41) for points 0(1) from the boundary: 

G c \ ose d(xi,yi;x + x 2 ,V2) = G closcd (0,0;a;, 0) 

_I i n (\ + X2 Z Xl \ - Vl ^ Vl + !) +2/2(^2 + 1) 

7T \ X J 2tt(x + X 2 — Xi) 2 

(xi, x 2 , j/i, and y 2 are all 0(1).) The Green function 
diverges as In 2, but if we are calculating correlations of 
LBMs along closed boundaries, we can use the arguments 
of section III to see that we only care about the parts 
of the Green function matrix that depend on both the 
row and the column indices. The part of Eq. (42) that 
depends on both (xi,yi) and (x 2 ,y 2 ) is 
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IX. THE ASM WITH DISSIPATION 



G c lo S ed{xi,yV, X + X2,V2) 



XlX 2 
TTX 2 



+ 



x\x 2 - x\x\ + x 1 y 2 (y 2 + 1) - x 2 y 1 (y 1 + 1) 



TTX° 



(43) 



Using logic identical to that in the bulk case, we can use 
the 0(l/x 2 ) part, to derive field identifications for LBMs 
along closed boundaries: 



- det(I + B„G MU ) (x T 1 B u x) 8686 (44) 

x is the vector of position coordinates parallel to the 
boundary. 

We now introduce a bond defect of strength fct> nd, 
along an open or closed boundary, between sites (q x i, q y i) 
and {q X 2 1 %2), that arc 0(1) apart. These bond defects 
can be analyzed as in the previous section. Along an 
open boundary, the bond defect is represented by 



-fcbondfei - q V 2) 8989 , 

7T 



(45) 



if q y i 7^ q V 2 — i.e. if the bond defect has a vertical com- 
ponent. On the other hand, if the bond defect is purely 
horizontal, it is represented by 



-fcbondfel + l) 2 (q x i - q x2 ) 2 d 2 9d 2 9 (46) 

7T 

A purely horizontal bond along an open boundary is rep- 
resented by a dimension four operator. 

For closed boundaries, the bond defect is represented 

by 



2 ~ 2 - 
khond{qxi - q X 2) 8989 , 

7T 



(47) 



if the bond has a horizontal component. If the bond is 
purely vertical, it is represented by 



^bond 

2tt 



(q yl (q yl + 1) - q y2 (q y2 + l)) 2 8 2 98 2 9 (48) 



Along closed boundaries, it is the purely vertical bonds 
that have dimension four. 

We have verified these field identifications of bond de- 
fects with more general calculations, involving multiple 
fields and multiple bond defects. This required gener- 
alizing Eq. (36) for multiple bond defects. However, the 
generalization is straightforward, and not particularly in- 
structive, so is not shown here. 



We now consider the addition of dissipation. As ex- 
plained in the introduction, this takes the ASM off the 
critical point, as shown by both numerical simulations, 
and an exact analysis [15-17]. The toppling matrix be- 
comes 



4 + t if i = j 
— 1 if i and j are nearest neighbors (49) 
otherwise 



Now, with each toppling in the bulk of the sandpile, t 
grains of sand are lost. 

If we set t = 0, we get back the original, critical ASM. 
It should be noted that the interpretation of the contin- 
uous t — > limit is potentially problematic. The mod- 
ification to A described in Eq. (49) is only sensible for 
integer t. It can be extended to rational t [16]. However, 
for t rational, but not an integer, the interpretation of 
the sandpile modifications (i.e. the B matrix) associated 
with a WAC is changed, so that what is meant by tak- 
ing a limit of infinitesimally small, rational t is unclear. 
However, presumably the fact that the massive results 
are good for all integers t > can justify an analytic 
continuation to t = 0. Regardless, we are certainly able 
to formally expand all correlation functions in Taylor se- 
ries about t = 0, which is what we do here. 

t defines an effective mass M for the sandpile, where 
t = a 2 M 2 [18]. a can be thought of as the lattice spacing. 
In looking at off-critical correlations of LBMs, we are 
interested in correlation functions where the number of 
lattice spacings between any two of the LBMs is 0{r/a). 
Taking the a — > limit then defines the way in which we 
simulataneously take t — > and distances between LBMs 
to infinity. 

For the off-critical sandpile, as discussed in [18], we can 
use the same methods as before to calculate correlations, 
with two modifications. First, we need to use a different 
B matrix than before. Previously, we required that our 
sandpile modifications be conservative, which meant that 
the each row and each column of B summed up to zero. 
However, for the massive sandpile we will often want to 
consider nonconservative B's. (For the unit height vari- 
able, the sandpile modification in Eq. (3) will no longer 
restrict the height of i to 1, but rather to any height from 
1 to 1 + 1.) We are most interested in LBMs associated 
with WACs. If we want to force the heights to the heights 
of the WAC, the B matrix must be changed to 



B = B c tB nc (50) 

B c is the B matrix that would be used for this WAC for 
the non-dissipative ASM (e.g. Eq. (3)), and 



1 if i = j , and i is in the WAC 
{B nc )ij = i height configuration (51) 
otherwise 
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By "in the WAC height configuration," we mean in the 
set of sites whose heights are fixed, and not in one of 
the bordering sites needed to to form the B matrix (e.g., 
for the unit height configuration, only one site is "in the 
WAC height configuration"). More generally, for other 
LBMs, B nc is the nonconservative part of the B matrix. 

Second, we need to use a new Green function between 
lattice sites. As always, the Green function is given by 
the inverse of the toppling matrix. In appendix A, we 
calculate the Green function in the limit t — > 0, when the 
distance between two sites scales as 1 j\ft, and find that it 
approaches (l/2ir)Ko(r), where Kq is the modified Bessel 
function of the second kind. The asymptotic expansion 
in \j\ft is then 



Go 



r cos r sin < 



Vt ' Vt 



2tt 



K (r) + Vfi(r, 



-t/ 2 (r,0)+0(i 3 / 2 ) (52) 



We have not calculated the / functions, because they 
turn out to not affect the universal parts of any corre- 
lation functions. In principle, they can contain bounded 
functions of 1/V, such as e tr ^, but we have not explic- 
itly indicated this i-dependence, since it does not affect 
our analysis. (Mahieu and Ruclle found, for <f> = and 
<f> = 7r/4, Eq. (52), and the specific forms of fi and / 2 [18]. 
They found that /i = for these angles, so it is possible 
that /i = for all </>, although we have not investigated 
this.) 

We decompose G uv as a sum of four N u x N v matrices: 



G„ 



'uv, J 



+ VtGuv.mw + VtG 



uv.co\ 



+ tG 



uu.both 



(53) 



G UVj j is a matrix in which every element is identical. 
G u „ ;row and G u „ ;CO i are matrices in which the elements 



depend only on the row index, or only on the column in- 
dex. Parts of G which cannot be written in these forms 
go into G^both- All four of these matrices are Taylor se- 
ries in V, whose 0(1) terms depends only on the Taylor 
expansion of (1/2tt)Kq (r), and whose higher order terms 
in \ft depend on /i, / 2 , etc. . . For example every el- 
ement of G„„ ;J is (l/2n)K (r uv ) + Vfi(r U v, (j> uv ) + . . .. 
The elements of G UVJOW and G uv ^ co \ depend on only the 
coordinates in the u th LBM, or on only the coordinates in 
the w th LBM, and thus require one derivative (finite dif- 
ference) of the Green function. The elements of G^both 
require two or more derivatives of the Green function. 

For correlations at the critical point, we saw that 
G mj j, G u „ jr0 w, and G UVyCO i could all be ignored in calcu- 
lating LBM correlations. However, the arguments there 
relied on the fact that every row and every column of 
every B u summed to zero. That no longer holds here, 
and we thus need to reconsider which terms of the Green 
function we need to keep. 

We expand the correlation functions in powers of i, 
and look for the lowest, nonzero, power of t. As with 
critical correlations, to prove the validity of Eq. (10), we 
need to show that the lowest-order term of the n-point 
function only comes from the parts of det(I+BG) with n 
terms off the block diagonal. However, the proof is much 
harder in this case; we sketch the proof in appendix B. 

Similarly to Eq. (7), we need to Taylor expand 



G [ -^=r„, 



COS( 



kl 'Vt ru 



. sm< 



f h -h 

(54) 

Then, just as in Eq. (12), we can write G uv as a sum 
of terms, each of which is the product of a length N u 
column vector and a length N v row vector. Defining \ u 
to be the vector of length N u , all of whose entries are 1, 
we have 



2ttG uv = K (r uv )l u l v 

+ [K' {r uv ) cos (puvVi^j {l u kl - fc„Lj) + (^K' (r uv )sm(t) uv Vtj - IvlC) 

- (K (r U v)t) ^sin 2 (f> uv k u k^ + cos 2 (j> uv l u l% + ^ sm(2<f) uv ) (kj% + /„fcj^ 

+ (K[l(r uv )t) (cos(2<^) (/"X - fc u fcj) - sin(20 ra ) (kj% + + • • • (55) 



In the ellipses, we have dropped not only terms of 0(t 3 ' 2 ) and higher, but all terms with /i or / 2 . The terms with f\ 
or / 2 are not necessarily higher-order in t than the terms shown, /i contributes terms of C(y/t) to lul^, and terms 
of 0(t) to l u ky ', ltJ„ , k u l^ , and l u l^ . / 2 contributes terms of 0(t) to l u lj. Similarly to Eq. (14), we represent G uv 
with a 3 x 3 matrix, , where each element of N„„ represents a different choice of row vector and column vector: 
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rv„, I., J.,, 



/ -t(sin 2 ^ o + cos(20)^')+O(^ /2 ) (t/2)sin(2^)(X -2^')) + O(i 3/2 ) -^cos0v^+) \ k u 
Kv = (i/2)sin(20)(A^-2A^) + 0(t 3/2 ) -i(cos 2 ^ o -cos(2(/.)^') + 0^ 3/2 ) -^0 sin M+O(*) f„ 

V ^cosM+0(t) i^sin<K/t + 0(t) X + O(t 1/2 ) / l u 

(56) 



To save space, we have here abbreviated 4> uv — > <fi and 
Ko(r uv ) — > i^o- Now, while some of the terms depending 
on /i and / 2 are ©(v 7 ^) and 0(t), they are higher or- 
der terms in G uv j, G UVtIOW , G u „ jCO i, and G uv ^oth, and 
thus are higher order in the specific matrix elements of 
N' uv that they contribute to. We will later see that this 
justifies dropping them. 

Just as in section III, when each off-diagonal Green 
function matrix is replaced by the product of a column 
vector and a row vector, each (I + B U G MU ) _1 B„ is brack- 
eted by a row vector to its left and a column vector 
to its right, producing a 1 x 1 matrix. We thus get 
a matrix, M' u , similar to the M u of Eq. (13). is 
3x3, rather than 2x2, because the vectors bracketing 
(I + T$ U G UU )~ 1 'B U to the left and to the right can now 
be k U7 l Ul or \ u . 

In principle, when calculating M^, B u and G uu should 
be the matrices for the massive sandpile. However, we 
only need most elements of in the limit t — > (this 
will be justified shortly). In this limit, we replace B u with 
B M C , and the elements of G uu with the normal, well- 
known, critical Green function. Thus, to lowest-order, 
the elements of M„ in arc unchanged: 



c 'u,kk = c uM + 0(t) 
c 'uM = + °(t) 

c'u,ll = c uM +0{t) 



(57) 
(58) 
(59) 



For the new entries of M' u , it is not hard to show that 



out to always be equal to the number of sites in the WAC 
height configuration [18]. We will see that our C is the 
same as their C, proving that their observation holds for 
all WACs. 
We now have 



Cu.kk + 0(t) c uM + 0(t) 
M' u = \ c uM +0(t) c uM +0(t) 
0(t) 0(t) 



The n-point correlation is given by 



0(t) 
0{t) 
tp u C + 0(t 2 ) 



(63) 



-Trace (M' N' M' N' M' N' ) 

— (((n— 1)! — 1) other trace terms) (64) 

It is easy to now verify that the higher-order terms in t in 
and N' uv that we dropped in Eq. (56) and Eq. (63) 
indeed give contributions of 0(t n+1 ^ 2 ) or higher to the 
n-point correlation, justifying the approximations used. 
We now have all correlations of LBMs in the ASM with 
dissipation. 

We can show that the correlations we have found are 
the same as the correlations of the field in Eq. (1). 
Mahieu and Ruelle proposed that the appropriate mas- 
sive extension of the LCFT in Eq. (20) is [18] 



M 2 

d z x{: 3039 : + — 



(65) 



This theory is still Gaussian, with correlation functions 



c 'u,lk — 


—irp 


11 


I + B u G Mtl 


J — 
c u,ll = 


-Pull 


I 


I + B u G Mtl 


c 'u,ll = 


-PuK 


I 


I + B u G Mtl 



B u k u - 0(t) 
BJ U = 0(t) 



- Pu ii{-m u , nc )i u + o{t z ) 

t Pu C + 0(t 2 ) , 



(60) 
(61) 



(62) 



where p u is the probability for the LBM, and C is defined 
as the number of sites in the WAC height configuration 
(as defined below Eq. (51)). c' u ll is the only element of 
M4 where we need the 0{t) term. Mahieu and Ruelle 
defined C as the coefficient in Eq. (1), and observed that 
for the 14 WACs that they considered, C always turned 



(6(z u )6(z v )) = K Q (M\z u 
(0(z u )6(z v )} = 
(6{z u )6{z v )) = 



(66) 
(67) 
(68) 



Other two-point correlations can then be obtained 
by taking holomorphic or antiholomorphic derivatives. 
Then, just as in section III, since the theory is Gaussian, 
we can write the n-point correlation of Eq. (1) exactly. 
The result is formally identical to Eq. (25), where the 
analogues of F and H in Eqs. (23-24) arc now 3x3 
(since the and 6 fields in Eq. (1) may have holomorphic 
derivatives, antiholomorphic derivatives, or no deriva- 
tives at all). Using the identifications in Eq. (26-28), and 
t = a 2 M 2 , the results agree with the correlation function 
in Eq. (64) (up to a proportionality factor, a 2n ). 
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APPENDIX A: MASSIVE GREEN FUNCTION 



where K is the modified Bessel function of the second 
kind. 

To find the leading-order critical limit, we take r — > 
(while rj\fl is still large), and use K Q (r) — > — ln(r), 
reproducing the first term of Eq. (6). 



The Green function for the ASM with dissipation is 
the inverse of the toppling matrix A, given in Eq. (49). 
In the limit of an infinite lattice, this can be found by 
Fourier transform [18]: 



G Q (m,n) 



v dpx r dpy 

_ w 27T J_„ 2tt : 



i(p x m+p y n) 



4 + t — 2 COSJXr — 2 cosp y 

(Al) 

We are interested in this integral either for m and n both 
0(1), or for m and n both proportional to 1/Vt, in the 
limit t -» 0. 

We saw in section IX, that for m and n both 0(1), we 
only needed G(m,n) — G(0,0), in the limit t — > 0. But 
this is just the standard, dissipation-free, lattice Green 
function, whose properties can be looked up in standard 
references — for example, [21]. So we only need to discuss 
the case where m and n are both proportional to \j\fl: 



t ^ I rcosft rsinc 
hmGo 1=- 

t^o \ y/t yft 



* dp x r d Py e ^Px cos 0+p B sin 0)r/v / t 



lim / ^± ^<L— (A2) 

t^oj_ 7[ 2ir J_ 7I 2ir 4 + t - 2cosp x - 2cosp y 

In the limit t — * 0, the exponential oscillates infinitely 
rapidly, so when multiplied by any function smooth in 
the limit t — > 0, gives an integral of zero. So no error 
is introduced by changing the region of integration to a 
small disc of radius e about (p x ,p y ) = (0,0). Similarly, 
because 



(A3) 



4 + t - 2 cos ^ - 2cosp y Pt+Py+t 



is smooth over this region, we can replace the first term 
with the second one. Changing to polar coordinates, the 
integral is 



lim 

t->o 



pe />2"7r 

/ Pdp 
Jo Jo 



pi(cos (ft cos a+sin (f> sin a)rp/ y/i 

da — „ (A4) 



(27r) 2 (p 2 + i) 
With a change of variables, the integral becomes 



,. „ / rcos0 rsin<; 
lim G 1=- 

t^O \ y/i VI 



_j_ r jdp_ r 2 * 

(2n)'J pi + lj 



(A5) 



APPENDIX B: THE TRACE FORMULA FOR 
MASSIVE CORRELATIONS 

In this appendix, we sketch the proof that the trace 
formula in Eq. (10) is valid for all n-point correlations of 
LBMs off the critical point. This requires showing that 
in the determinant dct(I + BG), if we only want terms 
up to 0(t n ), we never need to pick more than n terms 
off the block diagonal of I + BG. In other words, when 
we calculate the determinant with 



detX= Y, H P *i 



..X 



n,p(n) j 



(Bl) 



where p is summed over all permutations of 
{1, 2, . . . , \X\}, we never have more than n of the 
Xj jP (j)'s from the off-diagonal blocks. The off-diagonal 
blocks of 1 + BG all have the form B M G Mt) , u ^ v, where 
B„ and G uv can be written with Eq. (50) and Eq. (53). 
Since every row of B c u sums to zero, B c ,„G„„.j = 0, 
and B CiU G„„ ;CO i = 0. The off-diagonal block can thus be 
written 



-m nCtU G uvJ + 0(t 3 / 2 ) (B2) 

(We need the 0(t 3 / 2 ) terms for the calculation of the 
correlation functions in section IX, but do not need their 
explicit form for this proof.) Since there are terms of 
0(y/t) in the off-diagonal blocks, naively, to get the 0(t n ) 
contribution to the correlation function, we would need 
parts of the determinant with up to 2n terms off the 
block diagonal. So we need to explain why the terms with 
more than n terms off the block diagonal in fact have all 
their contributions to the 0(t n ) part of the correlation 
function cancel (as well as why all the terms with lower 
powers of t cancel) . 

We define a "row matrix" to be a matrix in which the 
entries depend only on the row index When we consider 
contributions to dct(I + BG), we consider, for each con- 
tributing matrix element off the block diagonal, whether 
it is from 1H cu G uv ^ OW i B c _ u G uv ,both? B nc ^G uv _j, or 
from the elements of 0(£ 3 / 2 ) or higher. We use several 
matrix theorems; the first two are general, not referring 
specifically to B or G. It is not hard to show 

Theorem 1 The determinant has zero contribution from 
terms that have both a matrix element from a row matrix 
in the (u,V\) block, and a matrix element from a row 
matrix in the (u, v 2 ) block. 
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By "zero contribution," we mean that while specific 
terms in Eq. (Bl) may be nonzero, when all such terms 
are considered, they cancel. Theorem 1 is easy to prove: 
if X a fj is the element from the (u, Vi) block, and X^s is 
the element from the (u, V2) block, we get a cancelling 
contribution from X a g and X 1 p. It is somewhat harder 
to show the following theorem, which we state without 
proof: 

Theorem 2 Suppose the (u,v) block (u 7^ v) is a row 

matrix, every column of the (v,v) block sums to one, and 
for every v! 7^ v, every column of the (v, u') block sums 
to zero. Then the determinant has zero contribution from 
the matrix elements of the (u, v) block; in other words, the 
determinant is unchanged if every element of the (u, v) 
block is set to zero. 

We note that B CiU G ut) row and B nc>u G ui;) j are both 
row matrices, and that every column of B CiU G ut ,^ row sums 
to zero, as does every column of B c U G U „ .both- Then, 
repeatedly applying theorems 1 and 2 allows us to prove 
the following: 

Theorem 3 In the determinant det(I + BG), suppose 
the (uq,U\) block (uq 7^ u\) has a contributing ma- 
trix element from B Ci „ G MoUl . row . Then, in the terms 
that produce a nonzero contribution to the determinant, 
there is an ordered sequence of distinct block indices, 
(ui, U2, ■ ■ ■ , u x ), x > 1, such that for all 1 < i < x, 
the (ui,Ui+i) block has a term from —iB nCtUi G UiUi+lt j. 
Furthermore, either 

1) The (u x ,u x ) block has a term of 0(t) or higher, 
or 

2) There is a term of order 0(t 3 / 2 ) or higher in an off- 
diagonal block, (u x ,u x+ i). 

A different (vq,v{) block (vq 7^ vi) with a matrix element 
from B C; „ G„ „ ljrow will produce an ordered sequence of 



distinct block indices, (vi,V2, . . .v y ), with no elements in 
common with (ui , 112 , . . . u x ) . 

Next, for any nonvanishing contribution to the deter- 
minant, we define 



ci = ff terms off the block diagonal (B3) 
c 2 = # terms off the block diagonal that 

are exactly C(t 1/2 ) (B4) 
C3 = # terms off the block diagonal that 

are C(t 3/2 ) or higher (B5) 
C4 = ff terms on the block diagonal that 

are 0(t) or higher (B6) 

Theorem 3 shows that each term in the determinant of 
type C2 can be associated with a distinct term of type C3 
or C4, so that C3 + C4 > C2. The number of powers of t 
from this contribution to the determinant is then at least 



(ci - c 2 - cs)(l) + 02(1/2) + c 3 (3/2) + c 4 (l) > ci (B7) 
So if we only want 0(t n ) contributions to the correlation 
function we should never have more than n terms off 
the block diagonal. Furthermore, we want a connected 
correlation function, so we should always have exactly n 
terms off the block diagonal. This concludes the proof 
that Eq. (10) is valid for off-critical n-point correlations. 

To get strict equality in Eq. (B7), we need c\ = n, 
C4 = 0, and C2 = C3. Furthermore, the terms of type 
C3 should be exactly proportional to i 3 / 2 . The fact that 
c 4 = means that in the diagonal blocks, I + B„G UU , we 
can set t = at the start of our calculations, as already 
seen by other means in section IX. 
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